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£ \ ABSTRACT 

o\ : 

CT) '. Neutron stars may be surrounded by fall-back disks formed from supernova 

^ ! core-collapse. If the disk circumscribes the light-cylinder, the neutron star will 

be an active radio pulsar spinning down under the propeller spin-down torque 

^3 ! applied by the disk as well as the usual magnetic dipole radiation torque. Evolu- 

tion across the P — P diagram is very rapid when pulsar spin-down is dominated 

q ! by the propeller torque. This explains the distribution of pulsars in the P — P 

! diagram. 
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•rH . 1. Introduction 

X' 

Efforts to understand the newly identified classes of neutron stars, in particular anoma- 
lous X-ray pulsars (AXPs - Mereghetti 1999) and soft gamma-ray repeaters (SGRs - Woods 
et al. 1999) have followed two avenues. Magnetar models, involving neutron star dipole mag- 
netic fields B ~ 10 14 — 10 15 G, above the quantum critical field B c = m 2 c 3 /eh = 4.4 x 10 13 
G, were advanced to explain the mechanism and energetics of soft gamma-ray repeaters 
(Thompson & Duncan 1995). Alternative models propose to explain the new classes of 
neutron stars in terms of conventional B ~ 10 12 G fields. These models involve accretion 
or propeller (Illarionov & Sunyaev 1975) torques from an accretion disk surrounding the 
isolated neutron stars (Alpar 1999, 2001; Chatterjee, Hernquist & Narayan 2000). Alpar 
(1999, 2001) argued that radio pulsars, dim thermal neutron stars (DTNs; Treves et al. 
2000), AXPs and radio quiet neutron stars (RQNSs; Chakrabarty et al. 2001) and perhaps 
SGRs represent alternative pathways of young neutron stars, distinguished by the history 
of mass inflow (M) from a fall-back accretion disk. This work classified young neutron 
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stars according to ranges of M, taken as representative constant values, with radio pulsars 
corresponding to zero or very weak M, such that the disk does not quench the pulsar magne- 
tosphere. The inner radius of the disk must lie at or beyond the light cylinder. Disks around 
radio pulsars were first proposed by Michel & Dessler (1981, 1983) and Michel (1988). For 
the AXPs Chatterjee, Hernquist & Narayan (2000) studied a fall-back disk with a specific 
time dependent mass-inflow rate M(t) taken to follow a self-similar thin disk solution, which 
entails a power-law decay of M(t) (Cannizzo, Lee & Goodman 1990, Mineshige, Nomoto 
& Shigeyama 1993). According to this model some neutron stars go through propeller and 
accretion (AXP) phases as M(t) evolves while others, at magnetic fields < 3 x 10 12 G become 
radio pulsars after a very brief initial accretion phase. 

What are the implications of a fall-back disk for radio pulsars? Marsden, Lingenfelter 
and Rothschild (2001a,b) argued that the presence of a propeller spin-down torque from a 
fall-back disk in conjunction with the magnetic dipole radiation torque may explain the large 
discrepancy (Gaensler & Frail 2000) between the real (kinematic) age and the "characteristic" 
age P/2P corresponding to pure dipole spin-down for the pulsar B 1757-24, and that pulsar- 
-SNR age discrepancies can be explained in a similar way. Menou, Perna and Hernquist 
(2001b) showed that this combination of dipole and propeller torques can explain the braking 
indices n < 3 of five young pulsars. 

In this letter we apply the combined spin-down torque, using the model of Menou, Perna 
and Hernquist (2001b), to the distribution of pulsars in the P — P diagram. Presence of 
torques other than dipole radiation torques, in particular torques that might arise from the 
ambient medium near a young pulsar in a supernova remnant was proposed first by Yusifov 
et al. (1995) to explain qualitatively the measured braking indices n < 3 of young pulsars, 
the discrepancy between real and characteristic ages and the distribution of pulsars in the 
P — P diagram. Gvaramadze (2001) invokes torques from circumstellar clumps in the SNR 
to derive a "true" age for PSR B1509-58 in agreement with the SNR age. In Section 2, we 
explore the combined torque model, using the canonical magnetic dipole torque together with 
the model of Menou, Perna and Hernquist (2001b) for the propeller spin-down torque, with 
a constant mass inflow rate M. Evolutionary tracks in the P — P diagram are presented in 
Section 3 along with analytical expressions for timespans of various phases, braking indices 
and other properties of the tracks. The P — P diagram is divided into strips delineated by 
tracks of constant magnetic field and mass inflow rate. Each strip is divided into period bins 
and a histogram is constructed for the number of pulsars in each period bin. A curve for the 
expected number of pulsars as a function of period is calculated according to the model and 
compared with the histogram for each strip. The results are discussed in Section 4. 
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2. Spin-Down of a Pulsar with a Fall-Back Disk 

We model spin-down of a pulsar under the combined action of magnetic dipole radiation 
and propeller spin-down torques as 

mn = -I3Q 4 - 7 (i) 

adopting the model of Menou, Perna and Hernquist (2001b) with their notation. The neutron 
star will continue to act as a radio pulsar as long as the fall-back disk does not protrude into 
the light cylinder. If the disk were detached from the light cylinder, it would not exert any 
torque on the neutron star and its magnetosphere. Assuming that the disk is attached to 
the light cylinder, the torque can be estimated as: 

N dlsk = -2Mrf c [n - Q K (ric)} = -2Mrf c H = -2Mc 2 /fi = - 7 /Q (2) 

where M is the mass inflow rate interacting with the light cylinder and being ejected from 
the disk; r ic = c/Q is the light cylinder radius, and rf c Q is the specific angular momentum 
extracted from the pulsar magnetosphere, since the Keplerian rotation rate VLk(j'Ic) in the 
disk is small compared to the rotation rate Q of the neutron star and its magnetosphere. 
The parameter 7 = 2Mc 2 = 2 x 10 31 M 10 erg/s is the rate of energy loss of the neutron star 
due to the propeller torque; M 10 is the mass inflow rate in units of 10 10 gm/s. The rate of 
energy loss due to magnetic dipole radiation is given by 

E dipole = -B ± 2 R 6 tf/Qc 3 = -(3Q\ (3) 

Here B± is the component of the dipole magnetic field at the neutron star surface in the 
direction perpendicular to the rotation axis and R is the neutron star radius. This defines 
P = 6.17 x 10 27 B ±:1 2 2 R 6 6 . Menou, Perna and Hernquist (2001b) give the solution of Eq.(l) 
for constant M as 

t = r[arctan((/?/ 7 ) 1/2 a 2 ) - arctan((/3/ 7 ) 1/2 ft(t) 2 ) ] (4) 
= r [ tt/2 - arctan((/3/ 7 ) 1/2 fi(t) 2 ) ] 

where Q iy the initial rotation rate of the pulsar, is always large enough to justify the second 
equation. The timescale r is r = I/2{jPf/ 2 = 4.5 x 10 7 M w 1/2 B^ 12 ~ l yrs. In the P - P 
diagram pulsars will follow tracks given by Equation (1), 

P = (4:7r 2 (3/I) P- 1 + (7/4tt 2 /) P 3 . (5) 

Menou, Perna and Hernquist (2001b) have used this model to explain the observed values 
of the braking index n = tlil/(l 2 . From their model we derived (M,B±) values (2.8 x 10 16 
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gm/s, 7.1 x 10 12 G) for the Crab pulsar, (1.5 x 10 16 gm/s, 8.9 x 10 12 G) for PSR B0540-69, 
(1.4 x 10 15 gm/s, 5.2 x 10 12 G) for the Vela pulsar, (2.8 x 10 13 gm/s, 8 x 10 13 G) for PSR 
J1119-6127, and (3.6 x 10 14 gm/s, 3 x 10 13 G) for PSR B1509-58. 1 



3. The Distribution of Pulsars in the P — P Diagram 

We used P and P data from the Princeton Pulsar Catalog (Taylor et al. 1996). Tracks 
corresponding to Eq.(5) are shown in Fig 1. At early times (small P), the dipole term 
dominates and the pulsar follows the left branch of the track, essentially spinning down 
under the magnetic dipole radiation torque. The braking index for pure dipole spin-down is 
n=3, and the left branch of each track follows the slope n' = 2 — n = — lin the log P — log P 
diagram, 

PP = (An 2 f3/I) oc B ± 2 . (6) 

The transition to the regime in which the propeller torque dominates occurs at the minimum 
of a track, at the period 

P = (27r)(/?/3 7 ) 1/4 = 27r/fio. (7) 

The slope n' = 0, and the braking index n = 2 at this point. The minimum value of the 
spin-down rate is found to be 

Pmin = (4tt 2 /3//)P - 1 + ( 7 /47r 2 /)P 3 = (8n / 1)^ /3f 4 . (8) 

The propeller spin-down torque is dominant along the right branch of each track. The 
asymptotic behaviour, reached already at periods of about 2P , has the slope n' — 3 (n — 
— 1). The evolution is very rapid on this branch. Neutron stars will spin-down to Q = 
(P = oo) in a finite time: 

t max = vrr/2 = 7 x 10 7 M^B^yrs. (9) 

A third of this time is spent until reaching the minimum P at P . About 0.23t max is spent 
until reaching n = 2.5, P = (3/7) 1 / 4 = 0.8P , when the pulsar has started deviating from 
dipole spin-down. The evolution slows down at periods around the turning point at Po to 
P = ^/3P , n = . From n = 2.5, P = 0.8P to P = ^/3P , n = it takes 0A4t max . As the 
pulsar proceeds into the propeller spin-down branch, the evolution speeds up, as P oc P 3 ; 
spin-down from P = y/3Po to P = oo takes only t max /3. 



1 The magnetic field of PSR Jl 119-6127 is B± = 8.2 x 10 13 G if one assumes dipole spin-down. Camilo 
et al (2000) report B± = 4.1 x 10 13 G in the discovery paper due to a missing factor of 1/2 in the definition 
of the neutron star's dipole moment. 
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Pulsar activity will turn off at reaching a critical voltage (the "death lines", "the death 
valley") at a period satisfying B_l j12 /P 2 — a — 0.1 — 1 (see Fig 1). The propeller spin-down 
tracks are parallel to the death lines. In our analysis of the P — P diagram we find that the 
minimum mass inflow rate for the pulsars in our sample is M min = 4.5 x 10 8 gm s _1 . The 
propeller spin-down track for M m i n bounds the pulsar population on the right and is of the 
form B^^/P 2 = 0.3. This track is already in the death valley. Pulsars with M < M m i n 
will reach the death valley and turn off while they are still on the dipole spin-down track. 
Pulsars with M > M min will evolve on their propeller spin-down tracks until they reach 
Pdeath = (B±, 12 /a) 1/2 at an age 

tdeath = t max [l - (2/tt) arctan(0.7a(M 10 r 1/2 )]. (10) 

This age for pulsar turnoff is close to t max for almost all cases. Even for the track at M min 
with a = 0.3 we find tdeath = 0.5t max - At times close to t max evolution is rapid. Few pulsars 
are observed at late times (long periods) along the propeller spin-down tracks. Therefore 
pulsar turnoff will not effect the distribution of pulsars in the P — P diagram except for the 
lowest M. The spindown lifetime t max is plausible in view of the low M values obtained for 
the pulsar tracks. A time dependent disk would reach M m 10 8 — 10 10 gm s _1 on timescales 
~ 10 8 yrs. If the disk is depleted before the pulsar completes its evolution the pulsar will 
switch to the pure dipole spin-down track corresponding to its magnetic field. 

The number AiV of pulsars with magnetic field B± and mass inflow rate M in a period 
interval AP at period P is 

AiV = (AiV(P;/3,7)/AP)AP = Rf(P,P)(At(P;{3, 1 )/AP)AP (11) 
= Rfi^P^PiP-Pn^AP 

Here R is the birthrate of pulsars in the galaxy, and f(P, P), the fraction of the galactic 
disk volume in which pulsars of period P and period derivative P are observable, describes 
selection effects. For simplicity we take f(P, P) to be constant, assuming that the variation 
in f(P, P) is negligible compared to the variation of P. We further assume that one and only 
one track (/3, 7) passes through each point in the P — P plane, such that the plane can be 
divided into strips separated by chosen (f3, 7) tracks. For each strip we count the pulsars in 
equal sized period bins and construct the histogram of AiV. We then choose a representative 
track ((3, 7) for the strip such that the function AiV(P) = constant /P(P;/3, 7) matches the 
histogram. 

The combined dipole and propeller spin-down model has a minimum P at the period 
Po on each track. AiV oc P _1 must have a maximum at P ; this is where a pulsar's motion 
across the P — P diagram is slowest. The model predicts that AiV oc P _1 will increase 
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in proportion to P at short periods, and decrease in proportion to P -3 at P > P . The 
histograms do have this property. For each strip, we fit the histogram by choosing ft to 
have the average value for the strip and choosing 7 such that Po is in the period bin with 
the maximum number AN max of pulsars. The normalization of the model curve is chosen 
to match the histogram at (Po, AN max ). Four strips separated by boundary tracks with 
-Bj_,i2 = 50, 4, 2, 0.8 and 0.1 are shown on the P — P diagram in Fig.l. For the upper four 
boundary tracks P = 0.3 s, and for the lowest boundary track at B ±:12 — 0.1, P — 0.8 s. 



4. Results and Discussion 

The number distribution of pulsars in the P — P plane is fitted remarkably well with 
the combined effect of dipole and propeller spin-down. The mass inflow rates range from a 
few 10 16 gm/s for the youngest pulsars to a minimum of 4.5 x 10 8 gm/s. The histograms 
for the four strips shown in Fig. 2 show reasonable agreement with representative model 
curves. The model curves representing the four strips are characterized by B± 12 = 35.5, 
P = 0.15 s (Fig.2a); B ±A2 = 3.16, P = 0.3 s (Fig.2b); B ± , 12 = 1.52, P = 0.3 s (Fig.2c) 
and P_i_,i2 = 0.57, P = 0.8 s (Fig. 2d). The agreement with the model is not sensitive to 
the choice of strip boundaries except for B^ 12 < 0.4 where a second peak appears in the 
histograms at P ~ 1.3 s, as seen in Fig. 2d. This excess may reflect disk depletion resulting 
in pulsars dropping down, from many propeller spin-down tracks at higher B^ and M, to 
their corresponding pure dipole spin-down tracks at periods close to death lines. 

Po is found to vary between 0.15 s and 0.8 s. From Eq.(7), the limited range of Po 
suggests a correlation between B and M. Such a correlation is not expected physically. It 
arises because a model with constant M is used to represent the real spindown which occurs 
under a time-dependent M(t). The period Po reflects the value of M at time to, when P 
is minimum. To understand the effect in a simple way, suppose at Po(to) the dipole and 
propeller torques experienced by the pulsar are equal. Suppose the mass inflow rate follows 
a power law decay M(t) oc t^ 11 . The pulsar is initially evolving down its dipole spindown 
track, with P(t) oc t 1 ^ 2 . We can define a Po(t) = P (/3,M(t)), the turning point into the 
propeller track, according to the current M(t). From Eq.(7), Po(t) = 27r( / 9/37(t)) 1//4 oc tf 1 ^. 
For early times P(t) < Po(t). One can imagine the pulsar moving along its dipole track 
with P(t) oc t 1 / 2 while the junction with the propeller track at the current M(t) moves 
down the dipole track, ahead of the pulsar, but at a slower rate Po(t) oc t^ 4 . At t — to, 
P(to) = Po(to) and the pulsar is just turning into its current propeller track. For t > t , 
the pulsar follows its current propeller track while the track itself continues to shift towards 
the bottom-right corner of the P — P diagram, as M(t) continues to decrease. To derive the 
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form of M(t) let us use the observation that Po(^o) = Po(P,M(t )) is similar for all pulsars. 
Taking P = 0.3 s gives M(t ) = 1.98 x 10 n S^ 12 . Dipole spindown (Eq.(6)) until t = t 
yields P = 1-23 x lO^E^^yr) 172 = 0.3 s or B 2 ± 12 = 5.95 x 10 6 /t (yr). Thus we find 
that M(t) = 1.18 x 10 18 (gm/s)t(yr) _1 , that is, rj = 1. For B ±jl2 = 1, t = 5.95 x 10 6 yr and 
M(t ) ^ 2 x 10 11 gm/s are obtained. 

The application here to all pulsars accross the P — P diagram and its success in fitting 
the pulsar distribution require and support the presence of low mass, low M disks, which are 
active for pulsar lifetimes of the order of 10 7 yrs and remain attached to the light cylinder. 
The time dependence of the disk and the mass inflow rate it supplies could be incorporated 
in a more detailed calculation. The self similar isolated thin disk models with a power law 
decay of M(t) were employed by Menou, Perna and Hernquist (2001b) who point out that 
for the Crab pulsar, PSR B0540-09 and PSR B1509-58, thin disks in the M ranges indicated 
by the braking indices of these pulsars are consistent with observational constraints in the 
optical though most of the disk luminosity would be in the UV where detections are much 
harder. For the Vela pulsar only a very small and highly inclined disk could be compatible 
with observational constraints. According to the classification of Alpar (2001) disks are 
present around all classes of young neutron stars, providing accretion for AXPs and acting 
as propellers on DTNs and RQNSs. For two AXPs and the RQNS in Cas A luminosities 
predicted by thin disk models are ruled out or tightly constrained by observations in the 
optical and IR (Coe & Pightling 1998, Hulleman et al 2000, Hulleman, van Kerkwijk & 
Kulkarni 2000, Kaplan, Kulkarni & Murray 2001). Further, it is likely that disks under 
propeller boundary conditions are not standard thin disks (Alpar 2001). As the mass inflow 
is stopped at the disk-magnetosphere boundary and largely ejected from the system, the disk 
may be enshrined in a corona and outflow of ejected matter, possibly with larger effective 
area and softer spectrum. Menou, Perna and Hernquist (2001a) have shown that the thin 
disks may become neutral and stop evolving as thin disks. The transition would quench 
the disk's luminosity, but it would also suppress the dynamical evolution of the disk if the 
viscosity of MHD origin is no longer operational. However, a smaller but nonzero viscosity 
should be operating in a neutral disk. Furthermore irradiation by the pulsar will probably 
keep the disk ionized and allow the mass inflow to continue, possibly as a power law decay. 
The minimum rotational energy loss rate of the pulsars is of the order of 10 30 erg s _1 . 
The luminosity of the irradiated disk may be shifted to different spectral bands (Perna and 
Hernquist 2000). 

The present application to the P — P diagram shows no indication of an effect on the 
pulsar distribution due to depletion of the disk before pulsar turnoff, except possibly for 
the oldest pulsars, well advanced in their evolution on propeller spin-down tracks, dropping 
down to pure dipole tracks near death lines. This may indicate that there is no separate disk 
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timescale. Disks interacting with pulsars via propeller torques may be driven by these torques 
to evolve on the propeller spin-down timescales of the pulsars. These problems, as well as 
the assumption that the disk remains attached to the light cylinder require further work on 
the coupling between the disk and the pulsar magnetosphere. A Monte-Carlo simulation of 
the P — P diagram will test the present model with independent random distributions of B 
and M and selection effects. 

We predict that some pulsars have braking indices < 2, like the Vela pulsar (though 
the n = 1.4 ± 0.2 braking index measurement for this pulsar should be treated with some 
caution as its timing behaviour is dominated by interglitch relaxation). About 2/3 of all 
pulsars will be on the propeller spin-down branch, with n < 2, at P > Pq. Half of these 
pulsars will have negative braking indices, — 1 < n < 0. Measurement of braking indices is 
unfortunately very difficult since the timing behaviour is dominated by noise (Baykal et al 
1999) and it seems also by interglitch recovery (Alpar and Baykal 2001). 

The model also predicts very small numbers AN oc P~ 3 of the oldest pulsars, which have 
long periods and large P values. These pulsars would give magnetic fields B± ~ 10 12 — 10 13 , 
perhaps 10 14 G with the combined dipole-propeller spin-down model. For these pulsars the 
pure dipole spin-down model would yield higher fields, extending into the magnetar range, 
and young ages. Thus kinematic age measurements of pulsars wih long periods and large P 
can distinguish between pure dipole spindown and propeller spindown. Propeller spin-down 
also increases the rate of energy dissipation in pulsars and in neutron stars evolving under 
propeller torques after pulsar activity is over. 
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Figure Captions 

Figure 1. The log P(s)— logP (s s _1 ) diagram. Strips used for the histograms in Fig.2 
are separated by the tracks shown, B ±jl2 = 50-4 (Fig. 2a), B ±:12 = 4-2 (Fig. 2b), B ±jl2 = 2 
- 0.8 (Fig.2c), and B ±jl2 = 0.8 - 0.1 (Fig.2d). The dashed lines are death lines B ±A2 /P 2 = a 
for a = 1 , 0.3 and 0.1 (left to right). 

Figure 2. Histograms showing the number of pulsars in period bins for each of the strips 
shown in Fig.l. The model curve on each histogram is drawn for the average (3 oc Bj_ 12 for 
the strip, with 7 oc M chosen such that the maximum of the histogram occurs at P = 
2vr(/3/3 7 ) 1 / 4 . 
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